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ABSTRACT 

To measure the stellar and orbital properties of the metal-poor RS CVn binary o Draconis (o Dra), 
we directly detect the companion using interferometric observations obtained with the Michigan 
InfraRed Combiner at Georgia State University’s Center for High Angular Resolution Astronomy 
(CHARA) Array. The P-band flux ratio between the primary and secondary stars is the highest 
confirmed flux ratio (370 ± 40) observed with long-baseline optical interferometry. These detections 
are combined with radial velocity data of both the primary and secondary stars, including new data 
obtained with the Tillingliast Reflector Echelle Spectrograph on the Tillinghast Reflector at the Fred 
Lawrence Whipple Observatory and the 2-m Tennessee State University Automated Spectroscopic 
Telescope at Fairborn Observatory. We determine an orbit from which we find model-independent 
masses and ages of the components (Ma = 1.35 ± 0.05 M 0 , Mb = 0.99 ± 0.02 M 0 , system age 
= 3.0 0.5 Gyr). An average of a 23-year light curve of o Dra from the Tennessee State Univer¬ 

sity Automated Photometric Telescope folded over the orbital period newly reveals eclipses and the 
quasi-sinusoidal signature of ellipsoidal variations. The modeled light curve for our system’s stellar and 
orbital parameters confirm these ellipsoidal variations due to the primary star partially filling its Roche 
lobe potential, suggesting most of the photometric variations are not due to stellar activity (starspots). 
Measuring gravity darkening from the average light curve gives a best-fit of /3 = 0.07 ± 0.03, a value 
consistent with conventional theory for convective envelope stars. The primary star also exhibits an 
anomalously short rotation period, which, when taken with other system parameters, suggests the 
star likely engulfed a low-mass companion that had recently spun-up the star. 

Subject headings: binaries: close - stars: activity - stars: imaging - stars: individual (o Draconis) - 
stars: variables: general 


1. INTRODUCTION 

RS Canum Venaticorum (RS CVn) stars are active bi¬ 
nary systems_with Ca H and K and photometric variabil¬ 
ity (jHalll fT976h . The components of these close binary 
systems are typically an evolved giant or subgiant pri- 
mary star with a s ubgiant or main s equence companion 
(lBerdvuginall2005t lSt,rassmeieill2009( f . 

In close binary systems, the stars experience changes 
in energy and angular momentum, circularizing the orbit 
due to forces from the aspherical mass distributions of the 
(partially) Roche-lobe-filling primary stars (IZahnlll9~77lf . 
Observations of open clusters have shown a transition pe¬ 
riod, the range of orbital periods between the longest cir¬ 
cular orbit and the shortest eccentric orbit, which mono- 
tonically increases with cluster age (Mazehi 120081 and 
references therein). Observations and predictions have 
shown that stars with orbital periods P 0 rb jS 10 days will 
have circularized while o n the main sequence, regard¬ 
less of spectral type (e.g.,lWalteHfl949l : lKoch fc Hrivnakl 
11981b iZahn fc Bouchetlll989fl . As stars evolve off of the 


main sequence and cool while becoming subgiants and 
giants, circularization is expected fo r stars with per iods 
Pnrh < 100 days (e.g. JMavor fc Mermilliodlll984 Iclaretl 
12009( 1. Therefore, RS CVn systems with orbital periods 
longer than 100 days are likely to retain primordial non¬ 
zero eccentricities. 

The single-lined RS CVn system o Draconis 
(omicron Draconis, o Dra, 47 Dra, HD 175306 , 
HIP 92512, HR 7125; G9TTT: IHerbig V SnaldindlU)^ 
Young fc Koniged 119771 : I Walter! 1 19851 iStrassmeier et al.1 
1989) has a previously-measured orbital period of 

Porb ~ 138 d ays with a n on-zero eccentric ity _ of 

e ~ 0.1 (e.g., Yound 119211 : iLucv fe Sweenevl 119711 : 
iMassarotti et al'l 2008(1 . The primary star of o Dra 
rotates faster than woul d be expected from tid al 
synchronization (e.g. IGIebocki fc S tawiko wskil 119881 : 
IMassarotti et alJl2008lf . 

To achieve this rapid rotation, the primary star must 
have increased its angular momentum. One possibil¬ 
ity is that the primary star engulfed a nearby com¬ 
panion, an event that would have had a strong im¬ 
pact on the primary star and its subsequent evolu- 
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tion. This evolution could shed light on the future of 
the Solar S ystem as the Sun expands to ingest plan¬ 
ets (e.g., ISchroder &; Connon Smith] I2008D . To in¬ 
vestigate this hypothesis, we aim to better dete rmine 
o Dra’s stellar parameters ( e.g., iBasri et all 119851 : 
iGurzadvan fc Cholakvanlll99 5j) and thereby resolve the 
binary system’s history. 

To advance our understanding of o Dra, we present 
the first detections of the low luminosity companion with 
six nights of interferometric data and the companion’s 
first radial velocity detections. We describe these data 
sets along with radial velocity and photometric data of 
the primary star in Section 2. We detail the analysis of 
observations with resultant orbital parameters in Section 

3. We describe the analysis of our photometry, show the 
detected ellipsoidal variations and eclipses, and measure 
gravity darkening of the primary component in Section 

4. We show our results on an Hertzsprung-Russell (H-R) 
diagram and discuss the system’s evolution in Section 5. 
Finally, we present the conclusions of our study of o Dra 
in Section 6. 

2. OBSERVATIONS 

2.1. Interferometry 

We obtained interferometric data at Georgia State Uni¬ 
versity’s Center for High-Angular Resolution Astronomy 
(CHARA) Array. The CHARA Array is a Y-shaped 
array of six 1-m class telescopes with non-redundant 
baselines varying from 34- to 331un l ocated at Mo unt 
Wi lson Observatory, California (Iten Brummelaar et al.1 
2005). Using all six telesc opes and the Michiga.n In- 
fraRed Combiner (MIRC; iMonnier et al.1 12004 120001 
IChe et a 1 .1 l2f)T~iTl . we obtained iJ-band (1.5 — 1.8 /im) 
data (eight channels across the photometric band with 
A/AA ~ 40) on UT 2012 May 9, 11, 12 and June 6, 8, 
17, 18; 2014 May 25, 26, 27, June 29, 30, and July 1. 

We detected the faint companion in the data from UT 
2012 May 9, 12; 2012 June 17, 18; and 2014 May 26, 27. 
The nights of observation without detections of the com¬ 
panion had insufficient uv coverage due to poor seeing, 
short observation lengths, or observations during eclipse 
leaving the companion undetected. We reduced and 
cali brated these data w i th the standard MIRC pipeline 

('seelMonnier et aT1l2007t[Zhao et ahll2009l IMonnier et al.1 

I2012L for pipeline details!. We used at least two calibra¬ 

tion stars for each night of data (see Table [[]). 

Three of our seven calibrators are A stars, which were 
revealed to be oblate due to rapid rotation during our 
analysis and required more information for calibration. 
Each of the stars were calibrated with a non-oblate star 
when possible, but some nights of observation required 
the use of other oblate stars. HD 185395 (9 Cyg) was 
used on 2012 Jun 17 to calibrate HD 192696 (33 Cyg) 
and HD 184006 (t Cyg). On nights when HD 185395 
was not available (2012 Jun 8 and 18), HD 192696 and 
HD 184006 were used to calibrate each other. HD 192696 
was used to calibrate HD 106591 (5 UMa) on 2012 May 
9, 11, and 12). The calibrated visibilities were fit to 
models of oblate stars to determine the mean uniform 
disk diameter, major-to-minor axis ratio, and position 
angle of the major axis (east of north), the mean values 
of which are in the footnotes of Table [0 


2.2. Radial Velocity 

We combined three radial velocity data sets of the pri¬ 
mary star and one set for the secondary to further con¬ 
strain our orbit of o Dra. 

Thr ee radial velocities published in [Massarotti et al.l 
(|2008f) were obtained with the CfA Digital Speedometer 
on the Wyeth Reflector at the Oak Ridge Observatory 
(Harvard, MA; 2004 September 1— 2005 April 14) . Two 
radial velocities also published in IMassarotti et ahl ((20081 ) 
were obtained with the CfA Digital Speedometer on the 
Tillinghast Reflector at the Fred L. Whipple Observa¬ 
tory (Mount Hopkins, AZ; 2006 June 14 — 2007 Septem¬ 
ber 21). Six additional radial velocities (2012 May 9 — 
2012 October 10) for the primary star were taken with 
the Tillinghast Reflector Echelle Spectrograph (TRES) 
at the Tillinghast Reflector (see Table [2]). 

With the Tennessee State University 2-m automatic 
spectroscopic telescope (AST) at Fairborn Observatory, 
AZ, we have determined radial velocities from 86 spec¬ 
trograms of the primary of o Dra taken between 2007 Oc¬ 
tober 11 — 2014 October 28 and 19 measurements of the 
companion radial velocity (2012 October 10 ^2014 Octo¬ 
ber 28; see Ta blcEH). See Appendix A. lFekel et al.1 (120091 ). 
and iRoettenbacher et al.l (120151) for details on these ob¬ 
servations. 

We additionally inc lude 18 a rchiv al radial velocities of 
the primary star from lYoungi (1 921 . 1902 July 15 — 1920 
August 9). 

2.3. Photometry 

o Dra has been monitored since 1992 with Tennessee 
State University’s T3 0.4-m Automatic Photometric 
Telescope (APT) at Fairborn Observatory. Our observa¬ 
tions span over 23 years from 1992 March 24 — 2015 May 
13, but with a gap during 2005—2011 (see Tableland 
Figure[l]). Our Johnson B and V measurements of o Dra 
were made differentially with respect to the comparison 
star HD 175511 (HIP 92594; B9; see Figure O and the 
check star HD 176408 (K1III). Details of the robotic tele¬ 
scopes and photometers, observing pr oced ure s, an d data 
analysis techniqu es can be found in iHenrvI (|1999t) and 
iFekel et all (120051) . 

2.4. High-Resolution Spectroscopy 

We use high-resolution spectra of o Dra covering the 
spectral range of 3700 — 7300 A obtained at the Nordic 
Optical Telescope (NOT) with the FIES high resolution 
echelle spectrograph. In the current work, the 1.3 arc- 
sec fiber giving a resolving power (A/AA) of 67000 was 
used. The exposure time for each spectrum was 180 s 
and resulted in a typical signal-to-noise ratio of 500 per 
resolution element (four pixels) at 6420 A. The observa¬ 
tions were carried out at thirteen epochs between 2012 
April 17 — August 15. All the spectra were reduced with 
the dedicated FIES reduction software FIEStool. 

3. ORBITAL ELEMENTS AND MASSES 

Using model fitting for the location of the unresolved 
secondary star with respect to the resolved primary star, 
we directly detected the companion of o Dra. Our models 
allow the primary star’s major axis and major-to-minor 
axis ratio, primary-to-secondary H-band flux ratio, and 
the secondary’s position to vary. During the fitting for 
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TABLE 1 

Calibrators for o Draconis 


Calibrator Name Calibrator Size (mas) a 

Source 

UT Date of Observation 

HD 106591 (5 UMa) 

oblate^ 

MIRC calibration 

2012 May 9, 12 

HD 125161 (t Boo) 

0.49 ± 0.03 

Bonneau et al. (2006) 

2012 May 9, 12 

HD 138852 

1.00 ±0.06 

Bonneau et al. (2006) 

2014 May 26 

HD 184006 (l Cyg) 

oblate c 

MIRC calibration 

2012 Jun 17, 18; 2014 May 26, 27 

HD 185264 

0.80 ±0.06 

Bonneau et al. (2006) 

2014 May 26, 27 

HD 185395 (9 Cyg) 

0.73 ± 0.02 

White et al. (2013) 

2012 Jun 17 

HD 192696 (33 Cyg) 

oblate d 

MIRC calibration 

2012 May 9, 12, Jun 17, 18 


Note. — a Some calibration stars are rapidly rotating and distorted from spherical. They were modeled as 
U-band uniform ellipses for calibration. For each star we list the iT-band uniform disk mean diameter (mas), 
major-to-minor axis ratio, and position angle of the major axis (°, east of north). 
b 0UD,mean = 0.78 =L 0.02, major/minor = 1.29 ± 0.05, PA ma j — 113 ± 5° 
c ^ud' mean = 0.72 ± 0.02, major/minor = 1.29 ± 0.04, PA maj = 92 ± 4° 
d ^UD.mean = 0.56 ± 0.02, major/minor = 1.41 ± 0.09, PA ma j = 115 ± 6° 


TABLE 2 

Radial Velocity Data of 
o Draconis (Wyeth & 
Tillinghast/CfA) 


HJD -2400000 

Primary (km s 1 ) 

53249.7015 

-27.27 

53306.5778 

-31.61 

53474.7086 

4.10 

53900.9303 

2.21 

54364.6628 

-31.30 

56056.9907 

-42.65 

56075.9832 

-30.30 

56110.9890 

3.86 

56132.6637 

-6.89 

56198.7292 

-41.88 

56210.5694 

-34.50 


Note. — Residuals on the pri¬ 
mary radial velocities are 0.39 km s — , 
The first five ve l ocities are published in 
iMassarotti et al.l (|2008D . 


TABLE 3 

Radial Velocity Data of o Draconis (AST/TSU) 


HJD -2400000 

Primary (km s 1 ) Secondary (km s *) 

54384.710 

-41.9 

56210.620 

-34.5 -1.2 

56211.589 

-33.4 

56212.589 

-32.6 

56227.713 

-12.2 


Note. — Table [3] is published in its entirety in the electronic 
edition. Residuals on the primary radial velocities are 0.22 km s —1 
and secondary are 1.6 km s — . 


the companion we weighted the closure phases ten times 
more strongly in the final y 2 than the squared visibilities 
and triple amplitudes in order to better identify asym¬ 
metries in the system to detect the faint companion (the 
detection of the companion is not sensitive to the factor 
of 10; see Figures [TU] — [13] in Appendix B). Errors for 
the primary star’s size and shape and the primary-to- 
secondary iJ-band flux ratio are based on the individual 
epochs. The positional errors of the location of the sec¬ 
ondary are error ellipses based upon the shape of the y 2 
surface used to detect the companion. 

The coordinates of the secondary detections on six 
nights (UT 2012 May 9, 12; June 17, 18; 2014 May 
26 and 27) are listed in Table [5] We measured the H- 



Fig. 1.— Johnson B and V differential magnitudes of o Dra 
acquired over 23 years from 1992 — 2015 with the T3 0.4-meter 
APT at Fairborn Observatory in southern Arizona. 



50000 52000 54000 56000 

JD-2400000 


Fig. 2.— Johnson B and V differential magnitudes of comparison 
star HD 175511 (with the check star 48 Dra). 

band uniform disk diameter of the primary star to be 
$ud = 2.115 ± 0.007 mas (limb-darkened disk $ld = 
2.189 ± 0.007, obtained by assuming a limb darkening 
power law exponent a = 0.27) with a major-to-minor 
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TABLE 4 

Differential Johnson B and V Photometric Data of o Draconis (APT/TSU) 


HJD—2400000 

A B (o Dra — 
HD 175511) 

AV (o Dra — 
HD 175511) 

A B (48 Dra“ - 
HD 175511) 

AV (48 Dra - 
HD 175511) 

48705.9454 

-1.119 

-2.330 

-0.101 

-1.282 

48706.9415 

-1.129 

-2.336 

-0.102 

-1.299 

48725.8916 

-1.136 

-2.350 



48732.8732 

-1.113 

-2.332 

-0.106 


48733.8831 

-1.118 

-2.334 

-0.105 

-1.287 


Note. — Table m is published in its entirety in the electronic edition. 
Q HD 175511 (HIP 92594), V = 6.94, B - V = -0.04. 
b 48 Dra (HD 176408, HIP 92997), V = 5.68, B - V = 1.16. 


axis ratio of 1.01 ± 0.03. The 77-band flux ratio for the 
primary star to the secondary is 370 ± 40, the highest 
confirmed binary flux ratio measured with long-baseline 
optical interferometry (RT Aur has an unconfirmed 77- 
band flux ratio of ~ 450:1, measured by MIRC at the 
CHARA Array; iGallenne et al.ll2015D . 

We calculate the orbital parameters of the binary by 
simultaneously fitting our interferometric and radial ve¬ 
locity data with Monte Carlo realizations. For our six 
interferometric points we present the error ellipses of the 
major and minor axis in Table [5] scaled to give our fit a 
reduced \ 2 = 1-00. The radial velocity errors are sim¬ 
ilarly scaled to require \ 2 = 1.00 (rmscfA.A = 0.39 km 
s , rms ast, A = 0.22 km s , rmsAST,B = 1-6 km s _1 , 
rmsYoun g ,A = 2.2 km s -1 ). 

The simultaneous radial velocity and astrometry 
Monte Carlo realizations of the orbit gave the orbital 
parameters and their 1 — a errors listed in Table [6] The 
orbit is represented in Figure [3] with the radial veloc¬ 
ity curve in _Figure [4j We use the conventions presented 
bv lHeint3 ( 1978 1. where the argument of periastron, u, 
follows the radial velocity orbit convention (the primary 
star with respect to the center of mass), which is differ¬ 
ent from the visual orbit convention (the secondary star 
with respect to the primary). The position angle of the 
ascending node (E of N), Q, is independent of definition 
being equivalent with respect to either the primary or 
secondary star. 

Our double-lined and visual orbit confirms previous 
analyses of orbital period and non-ze ro eccentricity (e.g., 
lYoundfl92H ILucv fe Sweenevlll971ll . while highlighting 
new evidence that the system is eclipsing (i = 89.6 ± 
0.3°). We determine the masses of the stars to be Ma = 
1.35 ± 0.05 Mq and Mb = 0.99 ± 0.02 M 0 ; the implica¬ 
tions to system evolution will be discussed in Section 5. 
Our orbit gives a par allax of 7r =_9.36±0.f0 mas, a value 
consistent with Ivan Leeuweiil (|2007l . 7r = 9.54 ±0.21 mas) 
but about 2 — o’ from the original Hi pparcos reduc¬ 
tion ( 7 r = 10.12 ± 0-43 m as; lESAIIl997h . confirmed by 
iPourbaix fc Boffinl (120031 7r = 10.27±0.42 mas). Adopt¬ 
ing the spectroscopic orbit presented here, the result¬ 
ing Hipparcos parallax becomes 7r = 10.15 ± 0.43 mas 
(Pourbaix, private communication). The origin of the 
2 — o’ discrepancy in parallax measurements is not well- 
understood; for our purposes in later sections, we proceed 
with our derived value of 7r = 9.36 ± 0.10 mas. 


4. LIGHT CURVE MODELS AND ELLIPSOIDAL 
VARIATIONS 



milliarcsecond 


Fig. 3.— Visual orbit of the RS CVn system o Dra with our 
dates of companion detection and their locations on the orbit (black 
error ellipses). The observed stellar radius of the primary star 
is plotted with a thick black line (o Dra A). The radius of the 
companion star, o Dra B, for the best-fit temperature of 6000 K 
is plotted as the small black circle. The light gray orbits represent 
fifty Monte Carlo realizations. Black lines connect the center of the 
detection error ellipse to the expected point in the best-fit orbit, 
which is overlaid in black (given in Table [6] with 1 — cr errors). At 
the southernmost point in the orbit, the secondary star is moving 
toward the observer. Note: the axis units are milliarcseconds (mas) 
with north upwards and east to the left. 


A previous study, IStrassmeier et alj (11989T) . suggested 
o Dra has photometric variations up to A V ~ O.Of mag. 
Other indications of activity include chro mospheri c 
Ca H and K emission and v a riation (lYoung fe Konigcs 
1977t ISi mon fe Fekell fl987l : IStrassmeier et al.1 1 1 9881) . 


Hall fe Persinged ( 198611 reported a possible rotation pe¬ 


riod of Prot = 54.6 days from photometric observa¬ 
tions. Using our NOT spectra, we can determine the 
value of usini. We modeled the Fe 6421 A line with a 
variety of v sin i models from ICastelli fe Kuruc3 (|2004T ) 
with Toff,A = 4530 K, logg A = 1.77, [Fe/H]= -0.5, 
microturbulence £ m = 1.4 km s , and macroturbu¬ 
lence £m = 6 km s -1 . The resultant value is usini = 
16 ± 0.5 km s -1 when combined with our observed pri- 















































Ellipsoidal Variations and Visual Orbit of o Draconis 


5 


TABLE 5 

Detections for the companion of o Draconis with respect to the primary 


UT Date 

JD —2400000 

Separation 

(mas) 

Position 

cr 

Error Ellipse 
Major Axis (mas ) 6 

Error Ellipse 
Minor Axis (mas ) 6 

Error Ellipse 
Position Angle (°) a 

Reduced 

x 2 

2012 May 09 

56056.97 

6.52 

23.1 

0.06 

0.04 

340 

1.3 

2012 May 12 

56059.95 

6.36 

23.7 

0.05 

0.03 

60 

1.5 

2012 June 17 

56095.87 

3.71 

203.1 

0.04 

0.03 

310 

1.6 

2012 June 18 

56096.82 

3.96 

202.3 

0.02 

0.02 

330 

1.5 

2014 May 26 

56803.87 

6.03 

202.5 

0.07 

0.05 

20 

0.9 

2014 May 27 

56804.85 

6.04 

202.0 

0.16 

0.04 

330 

1.1 


Note. — These detections give an //-band (1.5 — 1.8 /im) primary-to-secondary flux ratio for o Dra A to B of 370 ±40. The uniform disk 
fit for the primary star is 2.115 ± 0.007 mas (limb-darkened disk 2.189 ± 0.007) with a 1.01 ± 0.03 major-to-minor axis ratio. 

“East of North 

^Scaled error bars to ensure a reduced x 2 = 1-00 as described in Section 3. 


TABLE 6 

Orbital and Stellar Parameters of o Draconis 


Measured Parameters 

Value 

semi-major axis, a (mas) 

6.51 ±0.03 

eccentricity, e 

0.158 ±0.003 

inclination, i (°) 

89.6 ±0.3 

argument of periastron, to (°) a 

293.0 ±0.6 

ascending node, fi (°, E of N) 

22.9 ± 0.2 

period, P orb (days) 

138.444 ± 0.003 

time of periastron passage, T (HJD) 

2454983.0 ±0.2 

velocity semi-amplitude, Ka (km s — *) 

23.42 ± 0.05 

velocity semi-amplitude, Kb (km s — *) 

32.0 ±0.4 

system velocity, 7 (km s —1 ) 

-20.77 ±0.04 

H -band uniform disk diameter, #ud,A (mas) 

2.115 ±0.007 

H -band limb-darkened disk diameter, Old A (mas ) 6 

2.189 ±0.007 

primary major-to-minor axis ratio 

1.01 ±0.03 

B-band flux ratio, primary-to-secondary 

60 ±20 

V-band flux ratio, primary-to-secondary 

130 ± 80 

//-band flux ratio, primary-to-secondary 

370 ± 40 

rotational velocity, usin i (km s —1 ) 

16.0 ±0.5 

Derived Parameters 

orbital parallax, n (mas) 

9.36 ±0.10 

distance, d (pc) 

106.8 ± 1.1 

primary radius, Ra ( Rq) c 

25.1 ± 0.3 

primary luminosity, La (//q) 

220 ± 30 

primary surface gravity, log^A (cm/s 2 ) 

1.769±0.007 

primary mass, Ma (Mq) 

1.35 ±0.05 

primary rotation period, P ro t (days) d 

79 ± 3 

secondary radius, Rb (Rq) 

1.0 ± 0.1 

secondary luminosity, Lb ( Lq) 

1.3 ±0.2 

secondary surface gravity, logge (cm/s 2 ) 

4.43 ±0.09 

secondary mass, Mb (^ 0 ) 

0.99 ±0.02 

secondary temperature, T e ff b (K) 

6000±«o 

system age (Gyr) 

3.0 =F 0.5 

Literature Parameters 

primary effective temperature, T e ff a (K) e 

4430 ± 130 

primary metallicity, Fe/H/ 

-0.5 


Note. — “Radial velocity convention for primary with respect to the 
center of mass. 

b We applied a 3.5% correction from uniform to limb-darkened disk diameter, 
which is consistent with a limb-darkening power law exponent of a = 0.27. 
“Using limb-darkened disk diameter. 
d Assuming i vot = 2 or b- 

e T e ff, a is an average of tem pe rature s given by [C hristian fe Jane si 
(1977); Glgbocki &; Stawikowski (1977); Rutter & Schriiver (1987); 
Glgbocki & Stawikowski (1988); McWilliam (1990); | Luck| (|RI91|); 
Pourbaix & Boffin (2003); Bohp^VitensejJ|200jij) ; IMassarotti et al.| I2OO8I) ; 
Soubiran et al I 11201011 : IMcDonald et alj 1120121) . The 1 - a error is the 
standard deviation of these values. 

* [ Fe/H] is app r oxima t ed in stellar evolution models based upon value s given 
by IMcWilliaml dl990f) ; IMassarotti et al.| d2008D ; I Soubiran et all d2010|) . 
















































6 


Roettenbacher et al. 



Phase 

Fig. 4.— Radial velocity curves of the components of o Dra. 
The filled symbols represent measurements of t he prim ary star. 
T he fi lled diamonds represent observations from [Massarotti et al.l 
(2008)/CfA. The filled circles are the AST observations. The filled 
triangles are observations from lYound J1921J) . 1 — (7 errors in ve¬ 
locity are plotted unless the error is smaller than the symbols. 
Similarly, the open circles represent AST radial velocity data for 
the secondary star with 1 — cr error bars. The gray orbits are fifty 
Monte Carlo realizations (dark gray for the primary star and light 
gray for the secondary star) with the best-fit orbital parameters 
overlaid in black. See Table [ 6 ] for orbital parameters with 1 — cr 
errors. 

mary radius and orbital inclination (assuming z or b = ^rot) 
gives P rot ,a = 79 ± 3 days. 

We investigate the APT Johnson B and V differential 
light curves for evidence of starspots. We removed long¬ 
term variations (AU ~ 0.02) that may be attributed to 
axisymmetric spot structures or polar spot structures. 
We folded the adjusted light curves over the orbital pe¬ 
riod (P 0 rb = 138.444 days) and binned the data (0.01 in 
phase). The resultant Johnson B and V light curves are 
presented in Figure [5] The quasi-sinusoidal trend ob¬ 
served in the averaged light curves suggests ellipsoidal 
variations due to distortions of the primary star par¬ 
tially filling its Roche l obe p o tent ial (Pt, = 6 0.8 Pq, 
Ra/Rl = 0.42; e.g., [Roettenbacher et~aLll2015 i !l. Both 
light curves clearly show an eclipse (phase ~ 0.95) and 
evidence of a weak eclipse (phase ~ 0.40). Comparing 
the timing of the eclipses to the radial velocity curve and 
visual orbit, we see that the deeper eclipse is associated 
with the secondary star moving behind the primary star, 
revealing that the secondary star is hotter than the pri¬ 
mary star. 

The secondary eclipse provides flux ratios for Johnson 
B and V bands (60 ± 20 and 130 ± 80, respectively). 
Using the H- and P-band flux ratios, we are able to 
constrain the temperature of th e companion star. We 
use NextGen stellar atmospheres (lHauschildt et al.lll999f ) 
restricted to the bandpasses with the size and temper¬ 
ature of the primary star to determine the secondary 
star’s effective temperature (T e fl\B = 6000±3qq K), radius 
(Pb = 1.0 ± 0.1 M 0 ), luminosity (Lb = 1.3 ± 0.2 L 0 ), 
and surface gravity (logge = 4.43 ± 0.09 (cm/s 2 )). 

To model the observed ellipsoidal variations, we used 
the light-curv e-fitti ng software packag e Eclipsing Light 
Curve (ELC: lQrosz fc Hauschildtll2000l l. ELC accurately 



Fig. 5.— Differential folded and binned light curves of o Dra for B 
and V magnitudes plotted with error bars from the binning. Each 
data point is an average of data points spanning 0.01 in phase from 
the complete light curve folded on the orbital period. The quasi- 
sinusoidal signature of the averaged light curve is due to ellipsoidal 
variations caused by the primary star partially filling its Roche 
lobe potential. The lines represent the ELC models for ellipsoidal 
variations with the gravity darkening coefficient f3 = 0.00, 0.07, 
0.08, and 0.25, where /3 = 0.07 =b 0.03 is the best fit to the binned 
light curves. 

accounts for the star’s ellipsoidal shape which changes as 
the companion star moves along its eccentric orbit. We 
begin by modeling the light curves with no free param¬ 
eters. We assume that the orbital and rotational axes 
are aligned (z or b = Aot, -Prot = 79 days), although our 
results are not sensitive to the assu med P rn t- W e assume 
gravity darkening to be /? = 0.08 (jLucvl fl967h . The el¬ 
lipsoidal variations from the modeled light curves agree 
remarkably well with the observed light curve of o Dra 
confirming that the coherent quasi-sinusoidal signature 
is due to ellipsoidal variations (see Figure [5]). 

The long-term variations, ellipsoidal variations, and 
eclipses account for much of the large changes in the 
light curves, s uggesting that the previously identified 
stars pots fe.g.. lHall fc Persingerll 1 98fl iStrassmeier et al.1 
IT989h were instead observations of a combination of 
these effects. Additionally, the absorption lines of the 
NOT spectra do not show clear evidence of rotationally- 
modulated temperatures due to starspots (see Figure HO). 

However, weak spot signatures are occasionally still vi¬ 
sually present after removing the eclipse and ellipsoidal 
variation model from the observed light curve. To deter¬ 
mine a rotational period based on the spot signature, we 
perform a power spectrum analysis in which we sampled 
the light curve on a grid of one-day spacing (inserting 0.0 
on days without data). We found the strongest signature 
comes from a period of P ro t = 75 days (see Figure [7]) , 
slightly smaller than the estimated 79-day period based 
upon r;sini. While P ro t = 75 days is within 2 a from the 
spectroscopically determined rotation period, this small 
difference could be attributed to differe ntia l rot ation, of¬ 
ten seen with RS CVn stars (lStrassmeieill2009fl . 

Because the starspots are relatively weak, they do not 
strongly contaminate the phase-averaged light curve of 
o Dra, providing an opportunity for precision analysis. 
Using ELC, we measured the level of gravity darken¬ 
ing by modeling the system with no free parameters ex- 
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Fig. 6. — Portion of the o Dra spectrum containing two Fe I 
(6419 and 6421 A) lines. Eight spectra are over-plotted spanning 
JD 2012 April 17 — August 15, approximately two rotation peri¬ 
ods. Starspots moving across the stellar surface will manifest as 
features moving through some absorption lines. While the lines 
deviate from Gaussian profiles, throughout the rotation of o Dra, 
the absorption line cores do not vary in a periodic way, suggesting 
a featureless surface. 
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Fig. 7.— Power spectrum of the Johnson B light curve of o Dra. 
In the top panel, power spectrum of the B light curve with the 
long-term trend removed shows a peak at half the orbital period 
(69.2 days), consistent with the signature of ellipsoidal variations. 
In the bottom panel, the power spectrum of the B light curve with 
the long-term trend, eclipses, and ellipsoidal variations removed 
shows a peak at 75 days. 


cept th e gra vity darkening coefficient, (3 from T e g oc 
(Ivon Zeipell 119241) . We found that the best-fit gravity 
darkening coefficient for o Dra is (3 = 0.07 ± 0.03 (er¬ 
rors from bootstrapping ov er observational s e asons ), sim¬ 
ilar to recent fin ding s o f iDiurasevic et ahi (20061 (3 = 
0.06 ± 0.01) and iRoettenbacher et al.1 ( 20151 3 < 0.1). 
These measurement s are consis tent with the canonical 
value /3 ~ 0.08 from iLucvl (Il967h . but differs from an a l- 
ternative value /3 ~ 0.21 from lEspinosa Lara fe Rieutordl 

mm- 


5. HERTZSPRUNG-RUSSELL DIAGRAM AND 
EVOLUTIONARY HISTORY 


Fig. 8 .— H-R diagram for o Dra. The dashed and dot-dashed 
lines are the main sequence and post-main sequence evolutionary 
tracks for 1.35 M& and 0.99 Mq stars with [Fe/H] ~ —0.5, re¬ 
spectively d Potter et al.l 120081) . The gray regions represent our 
1 — <t mass errors with the solid black line representing the zero age 
main s eq uenc e. The dotted line is a 3 Gyr isochrone (PHOENIX; 
I Potter et al.I [20081) . The measured location of the primary with 
1 — <j errors is indicated by the plus sign. The companion loca¬ 
tion is indicated with the long-dashed line (with 1 — a errors in 
luminosity). 


To understand why o Dra has P ro t < Porb, we inves¬ 
tigate the evolution of the present stellar components. 
We plot the location of the components of o Dra on an 
H-R diagram (see Figure [5J using our measured stel¬ 
lar parameters determined from our orbit and flux ra¬ 
tios. We include the zero-age main sequence and Dart¬ 
mouth stellar evolutionary tracks (Fe/H = —0.5, q /Fe 
= 0.0, PHOENIX-based models. iDotter et all 1200 81/ for 
interpolated model masses (Ma,iu odei = 1-35 ± 0.05 Mq, 
Mb ,model — 0.99 ± 0.02 Mq). 

Our primary star detection falls on the 1.35 Mq 
evolutionary track with an estimated tempera¬ 
ture of 4430 ± 130 K (iChrist an fc Janesl 119771 
Glebock i fc Stawikowskil 119771: iRutter fc Schr iiverl 


1987t iGlebocki fc Stawikowskil 19881 McWilliaml 1199 


119911: iPourbaix fe Boffirt 20031: Bohm-Vitensel 


Massarotti et al.l 120081 : iSoubiran et al.l 120101 : 
McDonald et al.ll2012f) . The detection of the secondary 
with 1 — a errors cross the main sequence of the 
expected evolutionary track, but at the upper limit 
of our temperature range (P e ff,B = 6 OOO //300 K). Our 
system age estimate most strongly depends upon 
the mass of the primary star; for the primary mass 
{Ma ,model = 1-35 ± 0.05 Mq), we determine an age of 
the system of 3.0 =F 0.5 Gyr. 

We can use our knowledge of the evolutionary state of 
o Dra to investigate three possible explanations why the 
rotational period is faster than its orbital period. First, 
we conclude that the star could not have evolved right 
off of the main sequence. The rotational velocity is more 
rapid than expected (~ 3 km s -1 ) for the evolution of a 
1.35 Mq, early F main sequence star (based upon initial 
rotational pe riods of ~ 2 days and main-sequence radiu s 
of - 2 R q ; iNielsen et al.l 120111 iBovaiian et al.l Ft) 1 2 t) . 
Another component of the system must have imparted 
some angular momentum. 
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Second, secondary stars can spin-up primary stars 
faster than the orbital period in the case of eccentric 
orbits. Pseudosynclironous orbits are those that have 
the rotation of the star synchronized with the perias- 
tron passage of the secondary star (IHu 3 HMD. o Dra 
has a predicted pseud osync hronous rotation period of 
P ps = 123.3 ± 2.3 days (lHalllll986ll . which is much longer 
than our observed rotational period. Thus, the known 
companion of o Dra could not have spun up the system 
to its current rotational period. 

Third, a now-unseen companion could have spun-up 
the primary star. Because the secondary star has a mass 
of Mb = 0.99 ± 0.02 M 0 , the two objects that merged 
would have to have a significant difference in mass, other¬ 
wise the secondary star would be more evolved than the 
primary. With the primary and secondary stars falling on 
the same isochrone, the component that merged into the 
primary must have been only a small fraction of primary 
mass so as to not significantly affect its evolution. In or¬ 
der be consistent with rotational velocity observations, 
sufficient angular momentum would need to be imparted 
by the merging component to spin-up the ~ 1.04 Mq 
convective envelope of the 1.35 Mq primary staiQ- Un¬ 
fortunately, the rotation period of the star when the com¬ 
panion was engulfed is unknown and the primary star 
could have since slowed. While we cannot accurately 
estimate the mass of the consumed companion, we find 
from angular momentum arguments that the companion 
could range in mass from a giant planet to a low-mass 
star. If a low-mass companion were initially present in 
the system at 0.1 AU, it would be dynamically stable 
with the system’s stellar components on the timescale of 
the system age, 3 Gyr, but would be engulfed by the pri¬ 
mary star as the evolving stell ar radius approach ed the 
companion’s orbital radius (cf. iDavid et al.ll2003l l. 

6 . CONCLUSIONS 

In this work, we have made the first visual detections 
of the secondary star of o Dra using interferometric ob¬ 
servations. The iJ-band primary-to-secondary flux ratio 
is 370 ± 40, the highest confirmed flux ratio for a binary 
detected with long-baseline optical interferometry. With 
the astrometry and radial velocity data for both stars, 
we establish the first full three-dimensional orbit to de¬ 
termine orbital and stellar parameters. 

By folding and binning photometric data, we have 
shown evidence of ellipsoidal variations, gravitational 
distortions of the primary star caused by the close com¬ 
panion. The observed light curves are nearly identical to 
light curve models generated only from stellar and orbital 
parameters leading to the conclusion that the primary 
star has ellipsoidal variations, as opposed to long-lived 
starspots or active longitudes. After removing the model 
light curve, we observe only weak signs of rotationally- 
modulated starspots. Additionally, there could be poten¬ 


1 The mass of the convective en¬ 
velope was determined using EZ-Web, 

http://www.astro.wise.edu/$\sim$townsend/static.php?ref=ez 
K. H. D. Townsend’s Web-browser interface of the Evolve ZAMS 


tially active regions (e.g., axisymmetric spot structures 
and polar spots) that affect the stellar flux over longer 
periods of time with a global brightening or dimming. 
The folded, binned light curve also shows primary and 
secondary eclipses, which provide flux ratios to help con¬ 
strain the stellar parameters of the secondary star. 

Our new, high-precision orbital elements along with 
the folded light curve also allow for a measurement of 
gravity darkening. We find that /3 = 0.07 ± 0.03, a value 
of gravity d arkening consistent w ith conventional theory 
(lLucvlll967lf and previous results (|Diurasevic et a,l.ll2006l 
IRoettenbacher et al.ll2015f) . 

We have established that the primary star’s rapid ro¬ 
tation period could be due to the transfer of angular mo¬ 
mentum from a nearby companion. Specifically, a low- 
mass companion in a 0.1 AU orbit would impart sufficient 
angular momentum to spin-up the outer stellar layers be¬ 
fore being engulfed as the star ascended the giant branch 
while not dramatically altering the stellar evolution. 

In this series of papers, we have demonstrated that pre¬ 
cision interferometry at CHARA is capable of detecting 
the faint main-sequence companions of bright RS CVn 
primary stars. We are currently processing new data for 
another close, bright RS CVn system with the intention 
to publish these results in a follow-up paper. 
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Velocity (km/s) 

Fig. 9.— From a Fairborn Observatory spectrum of o Dra, the lower solid line is the average profile of the components summed over 168 
spectral regions. The upper line, arbitrarily vertically shifted, is the remainder after the average region around the primary component 
from all useful Fairborn spectra has been removed from the lower line. The position of the secondary is indicated by a tick mark below the 
summed residual spectrum. 


TABLE 7 

o Dra Orbital Parameters from Radial Velocity Curves 


Parameter 

Massarotti et al. (2008^ /CfA 

AST/TSU 

orbital period, P Q rb (days) 

138.45 ± 0.03 

138.46 ± 0.01 

center of mass velocity, 7 (km s —1 ) 

-20.8 ±0.3 

-20.87 ±0.02 

semi-amplitude, primary, Ka (km s —1 ) 

23.48 ±0.31 

23.37 ±0.03 

semi-amplitude, secondary, Kb (km s — ■*■) 


31.86 ±0.42 

eccentricity, e 

0.18 ± 0.02 

0.154 ±0.001 

time of periastron passage, T (HJD) 

2454982 ± 2 

2455812.8 ±0.2 

longitude of periastron, uja (°) 

292 ±5 

290.7 ±0.4 

Ma sin 3 i (Mg) 


1.34 ±0.04 

Mb sin 3 i (Mq) 


0.99 ±0.01 

a a sin i ( 10 6 km) 


43.97 ±0.05 

a,B sinz ( 10 6 km) 


59.92 ±0.79 

asinz (Rq) 

63.2 ± 0.9 


mass function, /(M) = (Mb sinz) 3 /(MA + Mb) 2 

0.177 ±0.007 


RMS velocity residuals, a a (km s — - 1 ) 

0.39 

0.16 

RMS velocity residuals, ctb (km s —x ) 


1.5 


APPENDIX 

A. AST RADIAL VELOCITIES 

The secondary lines of o Dr a were n ot o bviou s in the initial measurements of our individual AST spectra. Thus, 
similar to the method used bv lFekel et al.1 (|2015h . we subtracted the spectrum of the primary, which was obtained by 
averaging our spectra, appropriately shifted so that all the primary lines from spectrum to spectrum were aligned. 
This subtraction resulted in a very weak average summed profile of the lines in the residual spectra corresponding to 
the secondary component (Figure [9]). A Gaussian function was fitted to the very weak average secondary component 
to obtain its radial velocity. Because of the extreme weakness of the secondary lines and variable signal-to-noise ratio 
of the various spectra, most of the residual AST spectra did not produce measurable features of the secondary. 

For the primary we determined a usin i value of 15.6 ± 1.0 km s^ 1 , which is in good agreement with the value of 
16.0 ± 0.5 km s -1 from the NOT spectra. The equivalent width ratio of the average component lines, A/B, which 
corresponds to the continuum intensity ratio, was measured in several AST spectra that had the highest signal-to-noise 
ratios. We find a ratio of 128 ± 4 at a mean wavelength of about 6000 A, a value that is similar to the V band flux 
ratio of 130 found from the photometry. 

The parameters from the AST double-lined spectroscopic orbit and the combined lMassarotti et all (12008( 1/CfA orbits 
are presented in Table [3 
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Fig. 10.— Visibility curve of the 2012 Jun 18 observations of o Dra with CHARA/MIRC. The observed visibility curve is plotted in black 
with la error bars. The best-fit model from fitting for the companion is overplotted in as the black line (#ud — 2.106 mas, V(0) = 0.839, 
a = 0.00). The white line is our limb-darkened model (#ld — 2.189 mas, V(0) = 1.038, a = 0.27. The inset is the uu-coverage on the night 
of observation. 


B. INTERFEROMETRIC OBSERVABLES 

In Figures ITOl - [T3l we present a sample comparison of the calibrated o Dra data from 2012 Jun 18 and the best-fit 
detection of the companion from our x 2 -space fit. 
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Fig. 11. — Squared visibilitie s of the 2012 Jun 18 observations of o Dra with CHARA/MIRC. Each block of observations represents a temporal block of observations. The data and 
model are plotted as in Figure flOl 
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Fig. 12.— Closure phases of the 2012 Jun 18 observations of o Dra with CHARA/MIRC plotted as in Figure \TT\ 
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Fig. 13.- Triple amplitudes (multiplied by 10 4 for clarity) of the 2012 Jun 18 observations of o Dra with CHARA/MIRC plotted as in Figure ITT] 
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